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ABSTRACT 

The effects of anisotropic transverse collective flow on the HBT correlation function is 
studied. There exist three different physics contributions related to flow which affect the 
correlation function: anisotropic source shape, anisotropic space-momentum correlations 
in pion emission, and the effects related to the HBT measurement of the size of a moving 
source in different reference frames. Resolution of these contributions experimentally can 
lead to a detailed understanding of both collective flow in nucleus-nucleus collisions and 
the HBT technique itself. A method is presented which permits the derivation of model 
independent relations between the radius of a source measured in a frame in which it is 
moving and in its rest frame. 
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1 Introduction 



The discovery of anisotropic transverse flow in nucleus nucleus collisions at BNL AGS 
energies ]l| implies that some previous results should be reevaluated taking into account 
the effects of flow. One set of results in this category is the measurements of collision 
volumes using the Hanbury-Brown-Twiss (HBT) technique. In this paper we discuss 
how HBT results are affected by flow. In addition we show that the HBT study of 
nucleus nucleus collisions can provide valuable information on the understanding of 
flow itself. 

By anisotropic transverse flow we generally mean directed flow. We restrict our- 
selves to the consideration of symmetric collisions of identical nuclei. In this analysis 
we use the following geometrical definitions. In the transverse plane we define as 
the direction of transverse flow; "in - " is defined as the opposite direction; and u perp" 
is the direction perpendicular to the "in" direction. We use the coordinate system 
where the x axis coincides with and the z axis coincides with the beam di- 

rection. The reaction plane is then the plane defined by the "m" direction and the 
beam (the x-z plane). In this paper we consider only pion correlations (unless stated 
explicitly otherwise), and we call a pion source simply a source. 

We discuss the problem essentially on a qualitative level. Our goal is to find the 
physics which affects the HBT measurements, not to generate a complete set of formulae 
to describe the general case. In the discussion we keep in mind a heuristic picture of 
pion production and try to understand how different features of the production affect 
the HBT correlation function. In this very simple picture pions are produced from 
two sources different in origin. "Direct" pions are mostly produced in deep inelastic 
nucleon-nucleon collisions in the zone where the nuclei overlap. This source can be 
relatively small and is located close to the center of the line joining the centers of the 
colliding nuclei. The second source of pions in our picture is resonance (mostly delta) 
decays. This source size should be close to the freeze-out nucleon radius. In this naive 
picture it is clear that, for example, the "geometry" of the source of pions of different 
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rapidities (or different p t , or different orientation of pion momenta with respect to the 
reaction plane, etc.) can be very different. The effective source responsible for the 
emission of pions with a given rapidity can have nonzero longitudinal and transverse 
velocity. Note that we use the picture only to illustrate the sensitivity of the HBT 
function to different features of the source; one should remember that the real collision 
picture is much more complex. 

In general, the difference in the effective sources for pions in different rapidity (p t , 
etc.) regions can be treated mathematically as a correlation between the momentum 
and space position of pion emission. (Following many others we consider pion emission 
semi-classically, and do not discuss the space-time quantum mechanical uncertainty). 
For example, if low p t pions are produced mostly in the totally expanded stage of the 
collision, the extracted source size could be larger than that evaluated by using high 
p t pions, emitted from the hot stage. In Section ^| we discuss how these correlations 
appear in the expression for the correlation function. 

2 Geometrical shape of pion source 

We start with the question of whether it is possible to observe different source sizes 
looking at the emitting object from different directions. The answer is definitely "y es '\ 
for several reasons. One of them is simply the anisotropic source shape; due to trans- 
verse directed flow the effective source could be extended in the reaction plane. Directly 
produced pions are emitted mostly from the "center" (of the transverse plane), while 
pions from A (which undergo the collective motion to a greater extent) decays could 
have an "off-center" origin. Thus the source size appears to be different if measured 
from "in" or "perp" directions. 

We use RQMD 1.08 generated events to evaluate the approximate magnitude 
of the effect. We study production point distributions of pions with rapidity (in the 
laboratory frame) around yi a b ~ 3 which is close to beam rapidity in Au+Au collisions 
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for a projectile energy of 11.4 GeVV nucleoli. We select pions which are emitted in the 
transverse plane in the +x, —x, and y directions, and analyze their distributions in x, 
y, z, and t. Pions are considered as being emitted in +x direction if their transverse 
momentum components satisfy the condition < \p y \/p x < 0.5. All calculations are 
performed in the center of mass (Au+Au) frame. To suppress fluctuations related to 
very rare cases of decays of very long lived resonances, we consider only pions produced 
within the first 50 fm/c after the collision. The results are summarized in Table 1. 

From the first moments, one can see that the center of the observed pion source is 
shifted from the center line in the direction of flow. The variances of the distributions 
can be considered as squares of the effective source sizes in different directions. These 
are the sizes which one would measure if the pions were to carry the information on 
their production points (this is not the case in reality!). Due to the collision symmetry 
with respect to the reaction plane some coefficients in Table 1 (such as (y) for the +x 
and —x cases) are expected to be zero and others (such as (x) for +y and — y cases) 
are expected to be equal. The values actually obtained in the calculation can be used 
to estimate roughly the uncertainty in the results. 

The effective source velocities for each of the cases of +x, —x, +y, and — y directions 

of pion emission can be estimated from Table 1 using the formula: 

„ (rit)-(rj)(t) 

(t 2 )-(t) 2 ' {) 

which gives, for example, v\ +x ^ ~ 0.13, v\~ x ^ ~ —0.08, v^ +y ^ ~ 0.10, and for all 
direction of pion emission v z ~ 0.76. Note the difference between the magnitudes 
of y{ +x } and v\r x } which is just the difference in transverse flow velocities due to 
anisotropic flow. In collisions at this energy, flow is mostly carried by nucleons. Pions 
are involved in flow mainly through baryon resonance decays. It is important for our 
analysis that the pions from A decays almost do not "remember" the flow, since the 
change in pion momenta due to flow is small. The physical reason is the same as 
that for the low p t enhancement. When the resonance decays, the products share the 
resonance momentum in proportion to the ratio of their masses. The pion from A — > up 

4 



decay carries only about 1/7 of the momentum of the A. Thus we can expect for pions 
only about an extra 10 MeV/ c of transverse momentum in the direction of flow, if one 
takes reasonable values for baryon flow H]. Analysis of experimental data |7| yields an 
average value of v t ~0.3; for the longitudinal component it gives vi ong ~ 0.5. These 
values could be different from our estimate due to the particular region of pion phase 
space considered. 

The HBT analysis in principle permits to investigate the source size in many di- 
mensions, studying the dependence of the correlation function on different components 
of the relative pion pair momentum. Usually the so called "long-side-out" coordinate 
system is used for such an analysis. We remind the reader that the "long-side-out" 
coordinate system is defined in the following way (see Fig. [I]). "Long" is the direction 
of the beam; "out" is the direction of the pion pair transverse momentum; "side" is 
the direction in the transverse plane perpendicular to the "out" direction. In the case 
of anisotropic flow the picture becomes more complicated; one must also take into ac- 
count the direction of flow. Consequently, we denote, for example, by "m + -out" the 
"out" direction for the case when the pion pair momentum points in the direction of 
flow; and "m~-side" source size would mean the mean size of the source in the "side" 
direction as it appears for pions emitted in the direction opposite to flow. 

If one considers something like "m-out" and "m-side" source sizes, many interesting 
possibilities appear. One can see from Table 1 that, for example, the +x and — x 
effective sizes are quite different. One of the reasons for this is shadowing. Direct pions 
are produced most often in the region where colliding nuclei overlap. Let us look at the 
collision from the flow ("m + ") direction, where one expects the spectators and most 
of the nucleons of the projectile. Certainly these nucleons distort the real image of the 
pion source. From this direction one can see mostly the pions from the nucleon fireball 
and, consequently, one measures this fireball size. For an observer from the opposite (in 
the transverse plane) "m~" direction the source image is not distorted by shadowing 
so the picture could be quite different. From this direction one sees simultaneously the 
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hot pion fireball and, spatially separated from it, the nucleon fireball. The "out" size is 
the effective transverse size of the source in the direction of the transverse momentum 
of the pion pair; "m + -out" and "m"-out" sizes should be determined by the sizes of the 
two sources and their separation. This could explain the difference in "out" sizes from 
Table 1 seen from (V12.8 fm) and "m~" directions {yYJA fm). "Jn + -side" 

and "m~-side" sizes are the dimensions of the source perpendicular to the direction of 
the pion momentum as seen from and "m - " directions. The difference in "m + - 

side" and "m~-side" sizes (y/l9A fm and fm respectively) is mostly due to the 

difference in sizes of nucleon and pion fireballs, not to their separation. Another reason 
for the difference could be the shadowing of the pion source by nucleons in "in + " case, 
which effectively results in a larger observed source size. 

If this picture is correct one could study very interesting effects, varying the pseu- 
dorapidity of the correlated pair and the relative (azimuthal) angle between the pion 
pair and the reaction plane. In connection with the arguments made above it is also 
interesting to study in detail the pion triple differential distributions calculated with 
respect to the reaction plane.. For example, the shadowing discussed above could re- 
sult in nonzero values of the third harmonic Fourier coefficient of the pion azimuthal 
distribution. The rapidity (pseudorapidity) dependence, if observed, could provide 
information on the localization of dense nucleon matter. 

It is very interesting to compare the results of pion interferometry with the results of 
source size measurements by the two proton correlation technique. In our (oversimpli- 
fied) picture the proton source is different from the pion source. The effects discussed 
in the next two sections are also different for the two approaches, which makes the 
comparison more difficult but also more interesting. 
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3 HBT function and space-momentum correlations 

All of the arguments related to the difference between the and "in - " sizes given 

above can be formulated in terms of the correlation between the pion momentum and 
the space-time location of its production point. The reason for a separate discussion 
(in the previous section) is that the space-momentum correlations widely discussed 
in the literature are generally related only to isotropic transverse flow (any collective 
movement /phenomena in multiparticle production we call flow in this paper). The 
subject of our study is anisotropic flow. 

Below we rederive the expressions for effective source radii measured using the 
HBT technique, having in mind that this derivation can be helpful in understanding 
the flow contribution to the correlation function. Under a few simple assumptions the 
correlation function can be written in the "on-shell" form || : 



where S(x, p) is the source function, q = p\ — p 2 is the relative momentum, and 
P = pi + p 2 is the total momentum of the pair, P = (E,P); pi and p 2 are pion 
4-momenta. By consideration of the exact form of the correlation function one can 
generate the corrections || to the source parameters calculated with Eq. (fj). 

The mathematics of the space-momentum correlation is very simple; the correlation 
means that the source function does not factorize: 



A consequences of this is that the values presented in Table 1 depend on the pion 
momentum. We come to a trivial conclusion: the interferometry of pions with definite 
momentum is sensitive to the source which emits pions with just this momentum. The 
measured source sizes are not the sizes of the whole source but only of the effective 
region which emits the proper pions (the lengths of homogeneity ||). 

To give an idea of how the correlation function depends on the values discussed in 
the previous section and presented in Table 1, we derive below the expressions for the 




/ d A Xl d 4 x 2 S(x u P/2)S(x 2 , p/2) e -^' 1 - x2 ) 



(2) 




(3) 



7 



source radii. We consider the correlation function with a fixed value of P. In this case 
the correlation function depends only on three variables. From the fact that the pions 
are on the mass shell it follows that 



g = qP/E = qV, (4) 

where V is the velocity of the pion pair. Then, defining r = x\ — 22 = {t, r) the 
correlation function can be rewritten in the form: 

n( p s 1d /d 4 x 1 d 4 x 2 g(x 1 ,P/2)g(x 2 ,P/2)e- i ^ Vf - r ) _ 

C(q ' P) - 1+ LM%s(*,p/2)p • (5) 

The effective "mean square source size" can be defined through the second deriva- 
tive of the correlation function with respect to the component of q in the corresponding 
direction. There are diagonal: 

R 2 = -l/2(d 2 C( q )/d 2 q t )U =0 = m - r,) 2 )/2, (6) 

and cross terms [|J: 

4 = -l/2(9 2 C(q)/%^)| q= o = ((Vit - n){V 3 t - r 3 ))/2. (7) 
The correlation function in this case reads: 



C(q, P) = 1 + exp[- J2 &\ ~ 2 E 



ij) 



The interpretation of the expressions for the radii is straightforward: R 2 is the mean 
square of the distance between pions at the moment the second pion being produced 
(see also the discussion of this question in ||). 
In the "long-side-out" coordinate system 

q(Vt — r) = —q s ide r side ~ q<mt{Voutt — r out ) — qicmg^Ylongt ~ r long), (9) 

and it follows that 

2i&fa = (r 2 side ), (10) 
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2#L = ((rout ~ V out tf) = (r 2 out ) - 2(r out V out t) + ((V out t) 2 ), (11) 
2i?L s = ({non 9 ~ VW) 2 > = (rLg) - 2(r long V long t) + ((V long t) 2 ), (12) 

2Rlut,long = (( r long ~ V long t)(r out — V out t)) , (13) 

where we have introduced the notation 

J rf 4 xx d A x 2 f S(x u P/2)S(x 2 , P/2) 
U; [/d 4 x,S(x,P/2)] 2 ' 1 j 

In the absence of flow 

(r si de) = 0, (15) 

due to azimuthal symmetry of the collision ||. For the same reason R 2 side ^ ong = 
■Reside out = 0- For the case of non-zero flow these results are valid only for "in" sizes 
due to the symmetry of the collision with respect to the reaction plane. For the case 
of "perp" sizes all of the terms are in the general case nonzero. 

The expressions for the radii can be rewritten through the moments of one particle 
production point spatial distributions. For example: 

R 2 X = ((x-V x t) 2 }/2=(i(x 1 -x 2 )-V x (t 1 -t 2 )} 2 }/2 

= (( Xl - (x x )) 2 ) - 2V x ((x ltl ) - (x 1 )(t 1 )) + V 2 (( tl - ( tl }) 2 }, (16) 

and therefore the radii can be estimated using Table 1. The results of radii calculation 
are presented in Table 2 taking into account that in this case the mean transverse and 
longitudinal velocities of pions are approximately 0.38 and 0.87, respectively (in the 
center of mass of colliding nuclei). 

Due to the relatively small values of pion pair transverse velocities the radii R 2 and 
R 2 are not very different from the values of (x 2 ) — (x) 2 and (y 2 ) — (y) 2 in Table 1. They 
still resemble the features of the true source geometry. Due to the sizeable longitudinal 
pion velocity the values of R 2 are very different from the values of (z 2 ) — (z) 2 . Note also 
the possibility for non-zero values of cross-terms in the radii matrix. The relatively 
large value of R 2 XZ in the — x case could indicate spatial separation of the pion and 
nucleon fireballs. 
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4 HBT measurements of a moving source 



The values given in Table 2 are the parameters of the HBT correlation function as mea- 
sured in the center-of-mass frame of the colliding nuclei. As was shown in Section 2 the 
pion source in this frame has non-zero collective longitudinal and transverse velocities. 
Here we study how this motion distorts the HBT measurements of the source sizes. 

In some models it is possible to perform all necessary calculations in covariant 
form by introducing the 4-velocity of the source. However it is useful to derive model 
independent formulae by performing Lorentz transformations between different sys- 
tems. Our goal is to establish the relationship between the source size in its rest frame 
and the HBT correlation function measured in the frame where the source is moving. 
The correlation function by definition is the ratio of the two-particle invariant density 
and the product of the invariant one particle densities; it is invariant under Lorentz 
transformations. Using this property one can write: 

G(q,P)-G(q,P)-l + ifd*x>S(x>,P>/2)}i t 17 ) 

where the integrals are evaluated in the source rest frame, and the prime denotes the 
coordinate and momentum values in this system. 

Let us assume that our source moves with velocity v. We start with the case « 2 <1. 
In this case the momentum transformation equations are very simple: 

9o = <?o-vq = q(V-v), (18) 

q' = q- vg = q- v(qV), (19) 

P' = P - \E. (20) 

Then 

q'( x [ - x ' 2 ) = q ' x ' = ( q (V - v))t' - (q - v(qV))r. (21) 

One evaluates the second derivatives of the expression (pT7|) with respect to to obtain 
the radii. For example: 

R 2 X = -l/2(<9 2 C(q)/«9 2 ^)| q= o = ([t'(V x - v x ) -x> + K(rV)] 2 )'/2, (22) 
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R% = ([t'(V x - v x ) -x' + V x (r'v)][t'(V y - v y ) -y' + V y (rW)})>/2. 



(23) 



The other radii can be computed in an analogous way. The prime on the bracket (()') 
means that the mean value is evaluated using S(xi, P'/2) instead of S(x\, P/2). In fact 
the source function depends on P rather slowly and for an estimate of the radii one 
can neglect the difference for the case of small source velocity, which we consider here. 
Note that Eqs. 22 and 23 show explicitly the dependence of the correlation function 
radii parameters on source velocity v. 

The physical interpretation of Eqs. 22 and 23 is the same as for the case of the source 
being at rest, as can be shown by performing an approximate Lorentz transformation. 
The HBT correlation function measures the distance between the pions at the moment 
of production of the second pion, which is: 



as calculated to the first order of v. 

For the case of v <C 1 the values of the second moments of the spatial and temporal 
distributions are very close to each other in the both systems (note that in the source 
rest frame the quantities (tV) are zero). Taking this into account and using the 
values from Table 1 one can estimate the distortion of the correlation function due to 
transverse directed flow by considering Lorentz boosts only in the transverse plane. 
We find, for example, that about half of the difference between R 2 X in +x and — x cases 
is due to the transverse flow, with the remaining part attributable to the difference in 
source geometry for the two cases. 

We consider the case of arbitrary v and use the exact Lorentz transformations for 
one particular case when the flow velocity and pair velocity are directed along the x 
axis. In this case: 



tV-r = (lf + vr')V - (r' + vt') = t'(V - v) - r' + V(r'v), 



(24) 



qo - vq x q x V - vq x 



V-v 



(25) 



q 



q x - vq q x - vVq x 



= q x 



1-vV 




= q x 



(26) 
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where V is the pion pair velocity, and v is the flow velocity. From this formula and 
expressions for q' and q' x one can derive that: 

p2 , AV-v) x'jl-vV) ^, {{t'{V-v)-x'{l-vV)fy 

For a simpler case of instantaneous freeze-out (t = 0) and of a Gaussian source the 
analogous formula was derived earlier in [IIJ. Formula ( p^) is more general; it gives 
the dependence of HBT radii on the source velocity independent of any model. For 
anisotropic transverse flow the flow velocity depends on the orientation of pion mo- 
menta with respect to the reaction plane. This results in different apparent radii for 
the "m + ", "in - ", and "perp" cases even if the real source geometry is the same in all 
cases. 



5 Conclusion 

Anisotropic transverse flow produces anisotropy both in the pion source geometry and 
in the space-momentum correlation. It is also affects the interferometry size measure- 
ments because of nonzero effective velocity of the source. These three are the main 
phenomena responsible for the dependence of the HBT correlation function on the ori- 
entation of pion pair momentum with respect to the reaction plane. A detailed study 
of the experimental data and model predictions is necessary to disentangle all three 
effects, but such studies could be beneficial to our understanding of flow in nucleus- 
nucleus collisions and of the HBT technique itself. 
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Tables 



Table 1 : The first and second moments and correlation coefficients of spatial and temporal distribu- 
tion of charged pions in Au+Au collisions as seen from different directions. Pions are required to lie 
in the rapidity interval 2.7 < yi a b < 3.2 and have transverse momentum O.lAGeV < p t < 0.25 GeV. 
The impact parameter 3.0 fm < b < 6.0 fm. Units for the first moments are fm and for the second 
moments fm 2 . 
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Table 2: Matrix of source radii (in fm 2 ) found using the HBT correlation function measured from 
the different directions with respect to flow in the center of mass system of the colliding nuclei. The 
input parameters are taken from Table 1. 
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Figure Caption 

1. The definition of transverse momenta in the "long-side-out" coordinate system. 
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Figure 1: 
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